AD-A077  294 


UNCLASSIFIED 


MASSACHUSETTS  INST  OF  TECH  CAMBRIDGE  DEPT  OF  EARTH  A— ETC  F/G  8/5 
ANALYSIS  OF  LASER  RANGING  AND  VLBI  OBSERVATIONS  FOR  GEODETIC  PU— ETC(U) 
AUG  79  I  I  SHAPIRO  *  C  C  COUNSELMAN  F19628-77-C-0271 

AFGL-TR-79-0196  NL 


DDC  FILE  COPY  ADA077294 


A  PGI  ,  -TR -  79-0196 


ANALYSIS  OF  LASER  RANGING  AND  VLBI  OBSERVATIONS 
FOR  GEODETIC  PURPOSES 

1. 1.  Shapiro 
C.C.  Counselman 
T.  A.  Herring 


Department  of  Earth  and  Planetary  Sciences  — 
Department  of  Physics 
Massachusetts  Institute  of  Technology 
77  Massachusetts  Avenue 
Cambridge,  Massachusetts  02139 


( 


10  August  1979 


Final  Report 

1  October  1977  -  31  March  1979 


Approved  for  public  release;  distribution  unlimited 


This  work  was  sponsored  by  the  Air  Force  Geophysics  Laboratory 
and  Defense  Mapping  Agency  funded  under  contract  F19628-77-C-0271 


AIR  FORCE  GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
IIANSCOM  AFB,  MASSACHUSETTS  01731 


A 


Qualified  reau#>Q*/-> 

Tec  h  ni  c  al  ^fcTrrn  at  ion  *S  ervi  ce^  ®  ^  ®  ^  ^pply6 fQ  N  at  ionaf 608  6 


Unclassi f ied _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  hnlarad) 

J|-  j  1  /  REPORT  DOCUMENTATION  PAGE 


\2  GUV!'  ACCESSION  NO 


I  AFG^TR-7  9_2gl96J  j 

fl  TITLE  ltmU  Su  till  Hu) 

\  Analysis  of  Laser  Ranging  and  VLBI  _/> 

Observations  for  Geodetic  Purposes  .  j 


HEAD  INSTRUCTION: 

_ HKRORK  C'OMtM.KTING  l-ORM 

3  RECIPIENT'S  catalog  number 


f  s.  TYPE  OF  REPORT  A  PPWnTJ  rn/TRFD 

Final 

t  Oct  l-1977^|Mar  31,  1979 

A  PTLWE'U'WMIN^TJTi c/  H E pop T  NUW8FH  ’ 


V  AuTHORuJ 

I.  I. J  Shapiro^ 

C.  C . /Counselman 
I  T.  A.  /Herring 

"  'PTrerofi'M’N'G  organization  GameOTnd  address 

Massachusetts  Institute  of  Technology 
77  Massachusetts  Avenue 


’I  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Geophysics  Laboratory 
Iianscom  AFB,  MA  017  31 


V  NAM’  A  AuDRFsS'd/  tJitlereiU  from  iontrolUntf  Other > 


B  CONTRACT  ORTWlfNT  NUMHE 


F-^L  96  2 8 -7  7 -C- 027 1  v 

lOlOGBAMTLFMFUT. "pro  J|  :  1  P>' 

/  All,  A  A  l»«F>K  UNIT  N!J«Wh' 

f  /£  Jr  3  2  0  4pM AC  ^ 


IT,  61102F 


12.  REPORT  l)ATF 

0  August  1979 

13  N1IMAFB  or  PAGES 


IS  SECURITY  CLASS,  (ol  tiUs  report) 

Unclassified 


T  S/»  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


i6  distribution  statement  foi  this  Report) 

Approved  for  public  release;  distribution  unlimited 


I  1  7.  OlSTPtByTlON  STATFMFNT  ( of  t he  ehatract  "turret!  in  Work  JO,  It  <tlft>ront  from  Report) 


9  supplementary  notes 

This  work  was  sponsored  by  the  Air  Force  Geophysics  Laboratory 
and  Defense  Mapping  Agency  funded  under  contract  FI 9628-77-C- 

0271 

19  kFy  WORDS  (Continue  •«  rovorao  wide  1/  and  Identify  by  block  ntmtCMir) 

Very  Lone  Baseline  Interferometry,  Geodesy,  Intercontinental 
distance  measurement,  Polar  motion  and  Earth  Rotation 


3*.*  '• 


20  A  rE  ALT  '  Cnnttnu*  art  tayaraa  aula  II  tta >  an  rmry  mnd  Idantlfv  hy  hlorh  ntimhar) 

■^•From  threekVLBI  experiments  carried  out  in  1977  and  1978,  the 
distances  between  the  antenna  at  Onsala  (Sweden)  and  radio  tele¬ 
scopes  at  Haystack  (Massachusetts),  Green  Bank  (West  Virginia)  and 
Owens  Valley  (California)  have  been  determined  with  formal  stan¬ 
dard  errors  as  small  as  a  few  centimeters  and  repeatability  gener¬ 
ally  within  the  99?.  confidence  interval  of  the  root-sum-squaies  of 
these  standard  errors.  Overall  it  appears  that  suhdecimetei  pre- 
I  cision  has  been  achieved  in  these  estimates  ol  baseline  lengths. 


DD  ,V«m  1473  EDITION  OF  1  NOV  69  IS  OBSOL  f  TE 


Unc) assi t ied 


y  ^  /  /& 


SECURITY  CLASSIFICATION  OF  This  PAGE(»'Iiw  D»l«  Entered) 


20.  cont. 


Using  these  data,  we  have  also  determined  changes  in  UT1  and  in 
the  x  and  y  components  of  polar  motion. 
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I.  Coordinates  of  the  Antenna  i n  Onsala,  Sweden 

Using  the  data  gathered  in  three  VLBI  experiments  carried  out 
in  September  1977,  February  1978  and  May  1978,  we  determined  the 
coordinates  of  the  antenna  in  Onsala,  Sweden.  (See  Table  1  for 
experiment  details.) 

The  processing  of  the  1636  interferometric  group  delays  involved 
three  stages.  The  first  consisted  of  removing  group  delay  ambigu¬ 
ities  and  deleting  the  observations  for  which  there  was  no  signal 
detected  (based  on  a  small  correlation  coefficient  and  an  inconsis¬ 
tency  with  other  observations).  The  second  stage,  which  required 
a  large  amount  of  time,  involved  assessing  the  level  of  clock  and 
atmosphere  modelling  required  for  the  data  obtained  at  each  of  the 
observing  sites.  The  third  stage  was  the  generation  of  the  final 
solutions  in  which  we  obtained  estimates  of  baseline  vectors  and 
source  coordinates  using  individual  and  combined  data  sets. 

In  all  solutions  only  the  group  delays  were  used  to  estimate 
the  parameters.  The  delay  rates  were  used  as  a  means  to  assess  the 
quality  of  the  group  delay  solution.  This  procedure  was  found  to 
be  extremely  useful  in  deciding  on  the  level  of  clock  and  atmosphere 
modelling  required  for  a  given  data  set.  For  example,  if  the  addi¬ 
tion  of  extra  clock  and  atmosphere  parameters  improved  the  postfit 
residuals  of  the  group  delays,  while  degrading  the  delay  rate  resid¬ 
uals,  then  doubt  was  cast  on  the  reliability  of  this  solution. 

Table  2  gives  the  values  of  the  distances  between  the  inter¬ 
section  of  the  azimuth  and  elevation  axes  of  the  Onsala  antenna  and 
similar  reference  points  on  the  other  antennas.  The  results  pre- 


sented  give  values  from  individual  and  combined  solutions.  It 
should  be  noted  that  the  solutions  for  which  the  source  positions 
were  not  estimated  may  be  biased  due  to  the  greater  sensitivity  to 
source  position  errors  with  Onsala  included  in  the  experiments. 

The  self  consistency  of  these  results  is  gratifying  and  indi¬ 
cates  that  the  formal  standard  errors  quoted  may  not  be  much  smaller 
than  the  actual  errors.  The  solution  in  which  all  the  data  were 
combined  has  been  accepted  as  the  most  reliable  since  there  is  less 
possibility  of  clock  offsets  being  absorbed  into  the  2  coordinates 
of  the  antenna  positions. 

From  this  combined  solution,  our  estimates  of  the  baseline 
lengths  are: 

Haystack  -  Onsala  5,599,714.678  *  0.034  m 

NRAO  -  Onsala  6,319,317.765  +  0.037  m 

OVRO  -  Onsala  7,914,131.242  *  0.040  m 

These  values  are  believed  to  represent  the  first  intercontinen 

tal  distance  determinations  with  sub-decimeter  precision.  More  ex¬ 
periments  must  be  per formed, however ,  to  assess  the  repeatability  of 
the  determinations.  To  assess  accuracy  at  this  level,  comparisons 
would  have  to  be  made  with  determinations  based  on  an  independent 
and,  in  principle,  more  accurate  technique.  Unfortunately,  none 
now  exist.  The  best  such  technique  now  available  is  based  on  Dop¬ 
pler  observations  of  satellites.  Such  observations  have  already 
been  made  (V/.  E.  Carter,  private  communication)  ,  and  when  results 
become  available  for  these  baselines  comparisons  will  be  made. 


The  estimates  ot  tin  Onsala  antenna  coord  i  iki  t  es  (from  the  com¬ 
bined  solution),  for  the*  epoch  May  17,  J  •  >  7  K  a  <  . 

Onsala  X  3,  370,600.578  '  0.050  m 
Y  7  11,010.  ()!!'.  '  0 .  (I  !7  m 

Z  5 ,  34  0,837.  24  4  '  0 . 00  1  m 

where  the  coordinate  system  is  a  geocentric,  "rigid"  earth- I  i  :•<  :d 
system  with  the  Z  axis  parallel  to  the  mean  pole  of  rotation  ol 
1900  1905  as  defined  by  the  International  Latitude  Service  and  main¬ 
tained  by  the  Bill  (bureau  Internet  ional  (Jr  1  'lleure)  .  The  X  axis 
is  defined  to  bo  perpend  icular  to  the  axis  and  in  the  direction 
of  the  Greenwich  meridian.  The  Y  txis  completes  the  right-handed 
triad . 

Operationally,  the  location  ot  the  oriqin  ol  the  coordinate 
system  and  its  orientation  in  inertial  space  (mean  equinox  and 
equator  of  1950.0)  are  del  ined  by  the  coord  i  n.it  es  ol  the  Haystack 
antenna  reference  point,  which  wen  taken  to  be: 

Haystack  X  1,497,406.691  m 

Y  4,4  57,267.  130  m 

Z  -  4 , 296,882. 1 02  m 

ti  •  formulas  for  I  IT  I  ,  pole  positions,  sidereal  ti  me , 

nutation  and  precession  at  the  May  1978  epoch. 

In  earlier  reports*,  we  qave  values  ol  the  Onsal a-llaystack 
and  Onsala- NRAO  baseline  lenqths  determined  from  the  Si  ptember  1977 
lot  i  set.  When  the  earlier  determinations  arc'  corrected  for  a  change 

♦Quarterly  Status  reports  #1 ,  2,  and  3  lor  periods:  1  Oct  -  31  Dec 
1978,  1  .Jan  -  31  Mar  1978,  1  Apr  -  30  June  1978,  respectively. 


(i 


in  the  speed  of  light  from  299,792,5  km/sec  to  299,792.458  km/sec, 
the  differences  between  these  values  and  the  lenqths  given  in  this 
report  are  -0.5  and  -0.6  m.  These  differences  are  consistent  with 
the  quoted  formal  standard  errors  of  0.5  m.  The  earlier  report  also 
gave  the  coordinates  of  the  Onsala  antenna.  However,  a  comparison 
of  these  results  with  the  coordinates  given  in  this  report  is  not 
directly  possible  due  to  a  difference  in  the  (assumed)  coordinates 
of  the  Haystack  antenna. 

The  postfit  residuals  for  the  experiments  involving  the  North 
American  and  Onsiila  antennas  ranged  from  0.6  nsec  for  the  May  1978 
experiment  to  1.7  nsec  for  the  September  1977  experiment.  Figure 
1  qives  a  typical  postfit  residual  plot  for  the  May  experiment.  In 
all  cases  the  ratio  of  the  a  priori  estimate  of  the  accuracy  of  the 
observations  to  the  a  posteriori  estimate  was  of  the  order  of  1.2. 
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II.  Determination  of  Polar  Motion  and  UT1  Varia t i on s 

In  the  combined  solution,  where  all  observations  were  used 
simultaneously  to  estimate  baseline  vectors  and  source  coordinates, 
we  were  also  able  to  estimate  polar  motion  and  UTl  for  two  of  the 
experiment  dates.  The  third  experiment  date  was  used  as  a  reference 
since  VLBI  can  only  measure  changes  in  pole  position  and  UTl .  All 
four  stations  participated  in  the  May  1978  experiment  and  the  BIH 
values  were  taken  as  initial  values  at  this  epoch. 

Table  3  gives  the  differences  between  our  values  of  pole  position 
and  UTl  and  the  BIH  values  at  the  other  two  epochs.  Where  available, 
differences  deduced  from  other  techniques  have  also  been  included 
in  the  table. 

The  agreement  between  the  different  techniques  is  generally 
good,  although  in  some  cases  it  appears  that  the  formal  standard 
errors  may  be  smaller  than  the  actual  errors. 

The  largest  formal  standard  errors  associated  with  the  VLBI  de¬ 
termination  of  pole  position  are  in  the  y  component,  which  is  almost 
solely  determined  by  the  data  on  baselines  involving  Onsala.  This 
component  has  not  been  well  determined  because  we  have  only  limited 
data  on  these  baselines  in  the  September  and  February  experiments. 
However,  when  more  data  become  available  involving  Onsala,  we  can 
expect  to  have  the  y  component  of  pole  position  as  well  determined 
as  the  x  component. 

The  Lunar  Ranging  values  for  UTl  are  only  preliminary  and  were 
obtained  from  the  processing  of  2.5  years  of  data  from  November 
1975  to  May  1978.  The  average  formal  standard  error  of  the  results 
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was  '-0.5  milliseconds  (ms)  of  time.  However,  the  residuals  still 
show  systematic  trends  at  the  1  ms  level  which  are  thought  to  be 
due  to  deficiencies  in  the  mathematical  models  used  for  the  lunar 
orbit  and  rotation  (R.  W.  King  and  D.  S.  Robertson,  private  com¬ 
munication,  1979) . 


III.  Potential  Improvements  in  the  Dotermiruit ion  of  Source  Position 

VLBI  observations  involvinq  Onsala  and  the  North  American  an¬ 
tennas,  together,  potentially  yield  significant  improvements  in  the 
accuracy  of  source  positions  compared  to  the  accuracy  achievable 
with  the  North  American  antennas  alone.  This  conclusion  is  predomi¬ 
nantly  a  result  of  the  longest  baseline  length  being  doubled,  and 
its  declination  increased  from  6?7  to  11?0.  These  changes  afford 
the  improvement  in  resolution  shown  in  the  u-v  plane  diagram  of 
Figure  2.  Also,  the  '\,70u  difference  in  the  right  ascen¬ 

sion  of  the  Haystack- Onsaxa  and  the  Haystack  -OVRO  baselines  means 
we  are  able  to  obtain  the  right  ascension  anu  declination  of  a 
source  simultaneously.  Thus,  less  time  will  be  required  observing 
a  source  to  obtain  good  u-v  plane  coverage. 

Unfortunately,  the  number  of  observations  involving  Onsala  at 
this  time  is  still  too  limited  for  us  to  significantly  improve  the 
accuracy  of  the  source  positions.  This  condition  will  be  remedied 
when  observations  from  future  experiments  are  combined  with  the 
present  data  set. 


IV. 


Conclusion 


This  report  has  shown  that  subdecimeter  precision  can  be 
achieved  on  intercontinental  distance  measurements  using  VLBI.  It 
is  now  necessary  to  establish  if  we  have  achieved  the  same  level  of 
accuracy.  Comparison  with  satellite  Doppler  observations  is  one 
method  of  assessing  accuracy  and  it  is  hoped  that  this  comparison  will  be  possible 
soon.  A  partial  alternative  is  to  carry  out  further  VLBI  experi¬ 
ments  involving  Onsala.  In  November  of  this  year  a  Mark  III  experi¬ 
ment  using  Onsala  is  proposed.  This  should  provide  a  high  quality 
data  set  which  can  be  compared  with  the  May  1978  solution  to  obtain 
some  assessment  of  the  repeatibility  of  the  baseline-length  determi¬ 
nation.  This  does  not  provide  a  full  check  on  the  accuracy  since 
there  may  be  systematic  errors,  common  to  all  VLBI  experiments, 
which  are  currently  being  absorbed  into  clock  polynomials,  station 
coordinates  and  source  positions. 
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Table  1. 


Summary  of  VLSI  Experiments 


Involving 


Onsala 


a 


Dates 

21-25  Sept.  77 

24-26  Feb.  78 

17-19  May  197 

Experiment 

Duration  (hr) 

88 

50 

45 

Participating 

Onsala3  ^ 

Onsala3  ^ 

Onsala  ^ 

Stations 

Haystack0 

NRAOd 

Haystack 

OVROc 

Haystack0 

OVROc 

NRAOd 

Number  of  Group 

Delay  Observa¬ 
tions 

534 

272 

830 

Number  of  Sources 

Observed 

10 

9 

10 

Declination  Range 

-5.5  to 

2.3  to 

-5.5  to 

of  Sources  (deg) 

50.8 

50.8 

50.8 

66-ft  diameter  antenna  at  the  Onsala  Space  Observatory,  Onsala, 
Sweden . 

120-ft  diameter  antenna  in  Westford,  Massachusetts. 

130-ft  diameter  antenna  in  Owens  Valley,  California. 

140-ft  diameter  antenna  in  Green  Bank,  West  Virginia. 


e 


Onsala  only  participated  for  the  final  62  hr  of  the  experiment. 


is  less  susceptible  to  clock  offsets  being  absorbed  into  the  2  coordinate. 
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Table  3.  Pole  Posit  ion  and  UT1  Differences  Between  Estimates  from 
BIH,  VLBI  and  Other  Techniqu e s 


VLBI 

IPMS3 

Doppler*3 

LLRC 

24 

Feb.  1978 

x  component^ 

(m  arcsec) 

-33.6  ± 

2.2 

-19 

-60  ±  6 

_  _ 

d 

y  component 
(m  arcsec) 

15.1  ± 

7.0 

13 

-  6  ±  6 

— 

UTld 

(m  time  sec) 

-  3.0  ± 

0.1 

- 

- 

-3.2  ±  0.4 

24 

Sep  1978 

x  component*^ 

(m  arcsec) 

-80.3  ± 

7.7 

-62 

-53  ±  6 

- 

y  component^ 

(m  arcsec) 

-  4.9  ± 

9.1 

-18 

-20  ±  6 

— 

UTld 

(m  time  sec) 

-  1.0  + 

0.3 

-1.6  ±  0.4 

International  Polar  Motion  Service,  values  obtained  by  interpolation 
of  0.05  year  smoothed  values  compiled  by  D.  S.  Robertson  (NGS)  from 
the  Monthly  Notices  of  the  IPMS.  No  formal  standard  errors  are 
given . 

Naval  Weapons  Laboratory  and  Defense  Mapping  Agency  reductions  of 
Doppler  tracking  of  transit  satellites.  Values  interpolated  from 
0.05  year  smoothed  values  given  in  the  BIH  Annual  Report  for  1978, 
p.  D-13.  All  values  are  preliminary. 

Lunar  Laser  Ranging  preliminary  values  obtained  from  D.  S.  Robert¬ 
son  and  R.  W.  King  (private  communications) . 


d 


17  May  1978  taken  as  reference  point;  all  values  initialized  to 
agree  with  BIH. 


15 


5/18/78 
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The  letters  refer  to  observed  sources 

A=  NRAOl 50 ,  B=  3C84,  C*  4C3925,  D=  3C273B,  3C345  , 

F=  VR0422201,  G=  2134  +  00,  H=  3C454.3 

Lowercase  letters  denote  deleted  points.  (Achieved  by 
giving  the  observation  a  standard  error  of  1  sec.) 


Figure  2.  Comparison  of  L'-V  plane  covera 


